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We have recently proposed a novel self tuning mechanism to alleviate the famous cosmological 
constant problem, based on the general scalar tensor theory proposed by Horndeski. The self-tuning 
model ends up consisting of four geometric terms in the action, with each term containing a free 
potential function of the scalar field; the four together being labeled as the Fab-Four. In this paper 
we begin the important task of deriving the cosmology associated with the Fab-Four Lagrangian. 
Performing a phase plane analysis of the system we are able to obtain a number of fixed points for 
the system, with some remarkable new solutions emerging from the trade-off between the various 
potentials. As well as obtaining inflationary solutions we also find conventional radiation/matter- 
like solutions, but in regimes where the energy density is dominated by a cosmological constant, 
and where we do not have any explicit forms of radiation or matter. Stability conditions for matter 
solutions are obtained and we show how it is possible for there to exist an extended period of 'matter 
domination' opening up the possibility that we can generate cosmological structures, and recover a 
consistent cosmology even in the presence of a large cosmological constant. 

I. INTRODUCTION 

Over the last decade or so, as we have struggled to explain the nature of dark energy that is believed to be 
responsible for the observed acceleration of the Universe, interest has turned to the possibility that rather 
than being caused by an unknown form of energy density, the acceleration could be a result of a modification 
of Einstein's theory of General Relativity. It has resulted in an explosion of papers in the field, see [I] for 
a detailed review of the various approaches that have been adopted; one particularly interesting direction 
involves scalar-tensor combinations. It seems sensible to require that any theory maintains second order field 
equations, and the most general scalar-tensor theory satisfying that criteria was written down back in 1974 
by Horndeski [3] (it has recently been rediscovered in [5]). Such theories of modified gravity cover a wide 
range of models, ranging from Brans-Dicke gravity 6 to the recent models [7, 8J inspired by galileon theory 
[5]; the latter being examples of higher order scalar tensor Lagrangians with second order field equations. Of 
course all of these models can be considered as special cases of Horndeski's original action. Once the action 
was rediscovered, it did not take long before a perturbative analysis of the background evolution equations 
was carried out |10| lllj , which allows for a stability analysis to be performed on the various background 
solutions. 

In [5] we obtained a new class of models arising out of Horndeski's theory on FLRW backgrounds. The 
new models gave a viable self-tuning mechanism for solving the (old) cosmological constant problem, at least 
at the classical level, by completely screening the spacetime curvature from the net cosmological constant. 
In order to evade the famous no-go theorem of Weinberg [T^], the new solutions did not assume Poincare 
invariance to hold at the level of the solution (as Weinberg assumed), rather we allowed it to be broken in 
the scalar field sector. This is similar to a route adopted in [13] where the scalar field is allowed to break 
Poincare invariance on the self-tuning vacua, whilst maintaining a flat spacetime geometry. In [2] we provided 
a brief sketch of how the system works, showing that by demanding the self-tuning mechanism continues to 
work through phase transitions, so causing the vacuum energy to jump, we get powerful restrictions on the 
allowed form of Horndeski's original theory. Whereas the original model is complicated, with many arbitrary 
functions of both the scalar and its derivatives, we showed that by assuming matter is only minimally coupled 
to the metric (required to satisfy equivalence principle (EP) considerations) then once the model is passed 
through our self-tuning filter, it reduces in form to just four base Lagrangians each depending on an arbitrary 
function of the scalar only, coupled to a curvature term. We called these base Lagrangians the Fab-Four: 
Cj, C p , C g , C r , where the indices refer to John, Paul, George and Ringo. This was followed up in P3| with 
a detailed derivation of the conditions that lead to the four base Lagrangians just mentioned, in which we 
showed how they naturally lead to self-tuning solutions. Moreover in |14j we began to address the important 
question of the stability of the classical solutions to quantum corrections and demonstrated that at least 
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heuristically the self-tuning solutions can be guaranteed to receive only small quantum corrections thereby 
not spoiling the self-tuning nature of the solutions. 

The purpose of this paper is to begin the discussion of the cosmology associated with the Fab-Four 
Lagrangian. Without a sensible cosmology the model is nothing other than an interesting aside that may 
give some feeling as to how the cosmological constant can be addressed, but in itself does not have anything 
to say about our Universe. This is a non-trivial exercise. Note that some aspects of Fab-Four cosmology 
were touched upon in |15j . 

Recall from [5] that we are dealing with situations where the net cosmological constant may be large 
compared to any other energy density in the system. In the conventional cosmological scenario this would 
inevitably lead to a period of rapid acceleration, with no prospect of a radiation or matter dominated period, 
hence no chance for nucleosynthesis to take place or structures to form in our Universe. We will require the 
four potential functions in the Fab-Four action to act together and conspire to alleviate the influence of a 
net large cosmological constant before the final self tuning solution is reached. After all, we do not live in 
that solution just yet. To attack the problem, we will rewrite the equations of motion for the scalar field and 
the Friedmann equation in terms of a dynamical system allowing us to look for late time attractor solutions 
and to determine the stability of those solutions. 

The initial conclusions are positive, we are able to find combinations of the four potentials that do indeed 
lead to inflation/radiation/matter dominated like periods, with the latter two entering the self-tuning regime 
at very late times. For matter domination like behaviour, we explicitly show that there arc large regions 
of parameter space for which the corresponding solutions are perturbatively stable. These solutions are 
remarkable. Because we are interested in the case where the cosmological constant dominates the source, we 
have focussed on the case where no additional sources are present. Thus, we are able to find perturbatively 
stable matter-like solutions that are driven by a cosmological constant. What is happening is that the scalar 
field is working to screen the pressure component of the cosmological constant before its energy density. 
Eventually it will also screen the energy density, but the potentials allow for an intermediate period in which 
A essentially behaves like cold dark matter. This is the main result of this paper. 

The layout is as follows: in section [IT] we briefly recap the key Hamiltonian and scalar field equations of 
motion for the Fab-Four system arising from the original action of Horndeski [3] that is minimally coupled to 
matter. We do not rederive them, rather direct the reader to [14] for a rigorous derivation of the Lagrangian 
and evolution equations. We begin exploring the cosmology of this self-tuning scenario in section[TTlJ focussing 
on how each member of the Fab-Four behaves in isolation. To see how the various members behave in 
combination we rewrite the field equations as first order equations using a dynamical systems approach in 
section |IV| We switch off curvature in order to focus on the cosmological epoch prior to self-tuning, and 
find scaling solutions corresponding to different types of cosmology such as radiation domination, matter 
domination and inflation. Strictly speaking, some of these matter-like solutions are only fixed points for 
vanishing cosmological constant. Even so, as we show in section|Vj both analytically and through numerical 
simulations, they still provide an excellent approximation to the true cosmology even when there is a large 
non- vanishing cosmological constant. 

The reader not overly concerned with the details of how we arrived at interesting classes of Fab-Four 
potentials should probably skip section IV and proceed directly to section [V] Here we summarize the main 
findings of section [rv] as well as providing numerical simulations of solutions when spatial curvature is turned 
on. For the matter-like solutions, we also consider cosmological perturbations to weed out any problems with 
ghost and/or gradient instabilities. Whilst some solutions are unstable, others are perfectly well behaved. 
We conclude in section IvH 



II. THE SELF TUNING LAGRANGIAN - THE FAB-FOUR 



Given that we are interested in the cosmology associated with the Fab-Four, we will consider homogeneous 
and isotropic spatial geometries of the form, 



ds 2 



-dt 2 + a 2 (t) 



dr 2 
1 - kr 2 



\d0 2 



sm 



(1) 
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with k the constant denoting the spatial curvature. In |T3] we derived the sector of Horndeski's theory 
[3] that exhibits self-tuning. What does that mean? Given our cosmological background in vacuum ([lj, and 
the expectation that the matter sector can contribute a constant vacuum energy density, we identify with 
the net cosmological constant, p\ = P\ e + (pm)vac- In a self tuning scenario the net cosmological constant 
should not have an impact on the spacetime curvature, so whatever the value of p\, we still want to have a 
portion of flat spacetime 1 . What makes it applicable to cosmology is that this argument should hold when 
the matter sector goes through a phase-transition, changing the overall value of p^ by a constant amount, 
for example GUT phase transitions, EWK phase transitions etc. For this to work, any abrupt change in 
the matter sector has to be completely absorbed by the scalar field leaving the geometry unchanged. Hence 
the scalar field tunes itself to each change in p\ and this has to be allowed independently of the time (or 
epoque) of transition. The self tuning solution is Ricci flat, which tells us that at self tuning we have 

H 2 = -- 2 . (2) 

For k = we have a flat slicing of Minkowski, whilst for k < we have a Milne slicing. For k > no flat 
spacetime slicing is possible. With this demand of a viable self-tuning mechanism we were able to place 
powerful restrictions on the allowed form of Horndeski's original Lagrangian |14) . Whereas the original 
model is complicated, with many arbitrary functions of both the scalar and its derivatives, we showed that 
a self-tuning solution dramatically restricted the Lagrangian in form to just four base Lagrangians each 
depending on an arbitrary function of the scalar only, coupled to a curvature term. We called these base 
Lagrangians the Fab-Four, given by the following 



Cj = y/^gVMG^VvWvh ( 3 ) 

C P = V^rMP^VpWaWuVph (4) 

£ g = V^gVMR, (5) 

c r = V^gV r (4>)g, (6) 



where R is the Ricci scalar, G^ LV is the Einstein tensor, P^ap is the double dual of the Riemann tensor 
[I], and£ = W va & R^ af} - -iR^R^ + R 2 is the Gauss-Bonnet combination. Provided that {Vj,V p , V g } ^ 
{0, 0, constant} then these Lagrangian's naturally lead to the self-tuning solutions. Intriguingly, it follows 
that this constraint means that General Relativity is not a Fab-Four theory, which in itself is consistent with 
the fact that it does not have self-tuning solutions. By "self-tuning" , we mean that 

• the theory should admit a Minkowski vacuum for any value of the net cosmological constant 

• this should remain true before and after any phase transition where the cosmological constant jumps 
instantaneously by a finite amount. 

• the theory should permit a non-trivial cosmology 

Of course, this last condition ensures that Minkowski space is not the only cosmological solution available, 
which is important because we know from observation that the universe had to leave Minkowski space 
during its evolution. Fortunately this last condition is still allowed because the cosmological field equations 
are dynamical, with the Minkowski solution corresponding to a late time fixed point, meaning that once 
we are on a Minkowski solution, we stay there - otherwise we evolve to it dynamically. Note that for a 
homogeneous scalar, self-tuning is only possible for a Milne slicing of Minkowski. Indeed, the rate at which 
self-tuning kicks in and the solutions evolve towards the Milne Universe is controlled by magnitude of the 
spatial curvature, \k\. This will be evident from the numerical plots shown in section fVj 



1 Given the value of the tiny observed cosmological constant, for our purposes, flat spacetime is a very good approximation. 
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We begin our analysis of the cosmology by writing down the Hamiltonian density and scalar field equations 
of motion arising from the Fab-Four Lagrangian. Using the line element ([I]) the Lagrangians (|3| - ^ become 



Za 6 Vj 
-3a 3 K 



j,2. 



H 2 + k/a 2 ), 
'H{H 2 + k/a 2 ), 



C g = -6a 3 U 5 (0)(tf 2 -k/a 2 )-6a 3 V g ^{4>)H4>, 



Cr 



-ia 3 V r ^{cj))^H{H 2 + 3k/a 2 ), 



(7) 
(8) 
(9) 
(10) 



where "dot" corresponds to differentiation with respect to t, H = - is the Hubble parameter, and V r ^(4>) 



etc. Using 



+ 



dCj_ 

d<t> 



density in the presence of a matter source p m is 



T~L — T~Lj + Hp + T~Lg 



T~L r 



bare 

Pa — 



£j etc, we find that the Hamiltonian 



(11) 



where 



3H 2 



H L 



l~Lr 



-3K 



-6V g (<t>) 



H 5H 



H A + — 



-2W r , 4 (</>)<j>H [H 



H(f> 



MM 



v 9 {4>) 



The scalar equation of motion is given by 



E w = E 



E 



E 



E 



(4>) _ 



0. 



where 



E 



6 — 
dt 



3a 6 V H 



2 A, 



d 
dt 



a 6 V« 



■HA, 



3a 6 V„ 



'HA, 



-6- [^Vg^A,] +6a 3 V g ^)nSA 1 , 



-6a 3 V g ^)A 2 



= -24V^(0) 



dt 



and we have defined the quantity 



(12) 
(13) 
(14) 
(15) 

(16) 

(17) 
(18) 
(19) 
(20) 



A„ = H n - 



(21) 



which vanishes when we are on the self-tuning solution. As a result, it is easy to see that E^> also vanishes 
automatically during self-tuning. However, we note that the condition for self-tuning requires that the 
full scalar equation of motion should not be independent of a, and this is important as it ensures that the 
self-tuning solution can be evolved to dynamically, thereby allowing for a non-trivial cosmology. 

What is it we would like to recover from these equations? Ideally we want to find a cosmology consistent 
with observations, that does not rely on any particular value for the cosmological constant. It should be able 
to accommodate an early period of inflation driven by some combination of the four potentials, followed by 
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an extended period of radiation and matter domination during which nucleosynthesis could take place and 
in which structures could form. This would be followed by a late period of cosmic acceleration corresponding 
to the dark energy domination in which we find ourselves today. This is obviously a tall order, but as we will 
see, something that is not beyond the Fab-Four. In the next section, we will briefly examine the cosmological 
behaviour of each member of the Fab-Four in isolation to gain some intuition as to how each term will drive 



cosmology. This will be followed by a much more thorough analysis in section IV we introduce the powerful 
formalism of dynamical systems to rewrite the dynamics as a set of first order differential equations that we 
can then solve for their fixed points, allowing us to obtain a new set of cosmological solutions. 



III. THE COSMOLOGY OF EACH MEMBER OF THE FAB-FOUR 

To get a feel for how each member of the Fab-Four drives cosmology, we will briefly consider how they each 
behave in isolation, in the presence of a net cosmological, but no additional matter excitations. We neglect 
the latter because we are ultimately interested in the case where the net cosmological constant dominates 
the source completely, so any matter exciation will be subleading. As a result, we set p m = (p m )va,c in (H) 
and define the net cosomological constant pa = p^ re _+_ (Pm)\ 



It is convenient to rewrite the equations of motion ( II ) and ( 16 1 using N = In a as our evolution parameter, 

: d(f>/dN etc.), we find that 



as opposed to proper time. Using <j> = Hcj}' and a' = a, (where 



Hj + H p + rig + H r — —pA, 
Ef + + EW + - 0, 

where 



(22) 
(23) 



tin 



H r = -2AV r 



,2 H 2 3H 2 



13 tt4 



5H 



, H 2 H 2 



(24) 
(25) 
(26) 
(27) 



and 



E^HH 
Ef^H = 

eWaih = 



a 3 A 2 



3A 2 



a 3 / 2 (i?A 2 ) 1 /2 



H p a 9 / 2 (HA 2 y/ 2 A 2 
5A 2 -2^ 



-3^(A 2 a 
a 



4y 



(28) 

(29) 

(30) 
(31) 



Now, we are interested in the case where only one member of the Fab-Four is switched on. We therefore 
take Hi = —p\ = constant, and solve = for i = j,p,g,r. We know that at late times the solutions 
asymptote to a curvature dominated Milne universe H 2 — » — \ by the self-tuning mechanism [21 [14] . To 
examine what happens before that we consider the opposite regime in which the curvature is subdominant; 
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the results are listed below: 

'stiff fluid" (32) 
•'stiff fluid' (33) 
'radiation" (34) 
'curvature" (35) 

where M is some mass scale that arises as a constant on integration in each case. Vj and V p both behave 
like a stiff fluid with equation of state w = 1, which might have been expected since the corresponding 
Lagrangians contain derivative terms for 0, and we get a kinetically driven scalar field. In contrast, V g 
behaves like radiation whilst V r behaves like curvature. It is worth emphasizing that in each case the source 
is a cosmological constant with equation of state w = — 1, and yet the resulting cosmology behaves nothing 
like an inflationary de Sitter solution. The scalar field screens the pressure components of the source in 
a number of different ways, depending on which member of the Fab-Four is turned on. This allows us to 
be optimistic about extracting realistic cosmological solutions from the interplay between terms in the full 
Fab-Four theory. In the next section, we will show using a detailed dynamical systems analysis, that this is 
indeed the case 



Vj only : 
V p only : 
V g only : 
V r only : 



H z 
H 2 
H 2 
H 2 



M 2 

— (l + 0(k/a 2 )) 



a' 
M 2 

Ifi 
M 2 

If 
M 2 



(l + 0(fc/a 2 )) 



(l + 0(k/a 2 )) 



IV. A DYNAMICAL SYSTEMS APPROACH TO THE FAB-FOUR COSMOLOGY 



Our system of equations (11) and (16) are complicated second order equations and generally will not 
present tractable solutions. A powerful method to obtain attractor solutions is to adopt a phase plane 
analysis through a dynamical systems approach that allows us to both reduce the order of the differential 
equations but to also obtain the fixed point solutions without having to obtain the full dynamics of the 
system (see for example |16L I17j). This is a particularly powerful technique when we know what sort of 
fixed point solutions we are aiming for, in our case they correspond to radiation or matter domination or 
inflationary expansions. The reader who is simply interested in the cosmological solutions rather than their 
derivation may want to skip this section, and proceed directly to section |Vj 

We continue with N — In a as the "time" variable, and note that we will regularly come across the following 
combinations 

h= W a= l^' (36) 
from which it follows that when h = const we have 

H = H a h (37) 
a = a (t-t*)- 1/ft , (38) 

with i?Oj^o and t* as constants of integration; these will be used when we examine the fixed-points. If 
we were to neglect the spatial curvature then matter-like expansion would correspond to h = —3/2 and 
radiation-like expansion to h = —2. Curvature domination corresponds to h = — 1 and an inflationary 
trajectory would correspond to the region — 1 < h < 0, with the solution approaching de Sitter expansion as 
h — > 0. We will not consider the case where we have exactly h = 0, as it will correspond to a singular limit 
of our system. 

Another useful quantity when trying to understand particular solutions is the deceleration parameter that 
is given by 

q = 39 
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So for a ~ t p ~ t~ x l h we get 



Q = 



P(P ~ 1) 
p 2 



= -(fc + l) 



(40) 



The system of equations corresponding to (11) and ( |16| are quite complicated and so we require a few new 

(41) 



variables to fully establish the dynamical system. We do this by introducing 

x = H a 6' 



and 



Vp 



A, = 



V 



A„ = H~i* 



V p,<P 



V 



Vg 



Vr 



V 



g,<t> 



x r = m 



{V U f 



_ V P V P ,4>4> 



V a V 9 M 



,(j>Vr 

Z 1 '' = ~777 \2~> 



(42) 
(43) 

(44) 



where a, and 7^, (i = j,p, g, r) are constants. 

As in the previous section we specialize to the case where there is no matter excitation present, since we 
expect the net cosmological constant to dominate the source. Although the source has equation of state 
w = — 1, we will find that a judicious choice of Fab-Four potentials can mimic a standard cosmological 
evolution with any constant equation of state. In particular, we will show that even in the absence of an 
explicit matter fluid, and in the presence of a large cosmological constant p\, solutions exist that evolve as 
if the universe is matter dominated i.e. h = —3/2 in the language described above. 

Substituting (41 )-(44 1 into (11 1 and ( 16 ) we demand that the resulting system of equations are autonomous, 
which requires all the factors of the Hubble parameter H to scale out of the system. This follows from ( 11 ) 
where the term p\ is independent of H . Hence all the other terms need to be independent of H as well. 
Demanding that, we find the following relations 



7j =lp = lg 



P P 
Pg 

Pr 
- It 



4 - 2a, 
6 — 3a, 
2, 

4 - a, 
—a. 



The individual Hamiltonian parts are 



Hj = 3x 2 yj(3 — cr 2 ), 
H p = -3^(5 -3ct 2 ), 

Hn = ~ 



rtr 



-6y g [(1 - a 2 ) + X g x] , 
-24xy r (l-a 2 ), 



and combine in ( 11 ) to become 



H = Hj +H p + Hg+Hr + pa = 0. 



Similarly, the individual pieces of the scalar field equation of motion ( 16 1 are now 

E^/(3a 3 H a ) = 6xyj(l - a 2 ) + 2 (a - l)hxyj(l - a 2 ) + 2[xyj(l - a 2 )}' - X 3 x 2 yj {l - a 2 ), 



E^/(3a 3 H a ) 
E^/(3a 3 H a ) 

E^/(3a 3 H a ) 



= -9x 2 y p (l - a 2 ) - 3(a - \)hx 2 y p {l - a 2 ) - 3[x 2 y p (l - a 2 )}' + X p x 3 y p (l 
= -6X g y g (l -a)+ 2(1 - a)hX gVg (l -a)- 2[X g y g (l - a)}' 
+2p g X 2 g y g x(l -a) + 2X g y g (l - a) 2 , 



-My r (l + h) 



-a 2 (l 



8y r 



-a 2 {l-a) + -(l 



(45) 
(46) 
(47) 
(48) 
(49) 

(50) 
(51) 
(52) 
(53) 

(54) 

(55) 
(56) 

(57) 
(58) 
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and give the full scalar equation of motion 



E W = E W + E W + E W + E W = 0. 



(59) 



For completeness we can write down the individual pieces of the Euler-Lagrange equations for the scale 
factor, Ej^ = *^ — (l^r) W e * c ' wri ich combine to give the full scale factor equation of 



motion 



E {a) = E 



(a) 



E (a) + E (a) + E (a) _ 3a 2 flA 



0. 



(60) 



em 

Ei a) 



3a 2 {-4x 2 yj + 2hx 2 yj - 2(x 2 yj)' + x 2 y J (l - a 2 )} , 
'da 2 {2x 3 y p (3 - a 2 ) - hx 3 y p (3 - a 2 ) + [x 3 y p (3 - a 2 )]'} , 
3a 2 {4y' g + 4y g (2 - h) + 2{\ g xy g )' - 2h\ g xy g - 2y g (l + a 2 )} , 
3a 2 {I6xy r (l - a 2 ) - 8hxy r {l - a 2 ) + 8[xy r (l - a 2 )}'} . 



(61) 
(62) 
(63) 
(64) 



Of course, although we have written down the equations for the Hamiltonian constraint, the scalar field and 
the scale factor, we only require two of these sets as the third can always be obtained as a combination of 
the other two. 

We have nearly completed our derivation of the dynamical system. All that remains is to consider the 



evolution of the y, A and fj, coefficients in (41 1 - (44) and a in (36). These are obtained directly from their 
definitions. All the A's evolve as 



A- = [-ah + Xix(fj,i - 1)] A,, i = j,p,g,r. 

The y's evolve according to 

y'i = \Pih + Xtx) yi, i = j, p, g, r 

whereas a evolves according to 

a' = -{l + h)<7. 

We note that fixed point solutions for Aj and yi correspond to 

ah — Xix(fii — 1) 
Pih = —XiX 



(65) 
(66) 
(67) 



(68) 
(69) 



Similarly from (67) fixed point solutions for a exist for h = — 1 and a = 0. Our complete system of equations 



for investigatating the cosmology is therefore given by (54), (60), (59), (65), (66) and (67). We also have the 
definitions d36b , l4l| - fl44l and the conditions J45l - M9T 



In order to obtain closed form solutions, it is convenient to assume that all the pbi are constant, which 
generically corresponds to power law potentials 



Vi - 



j,p,g, 



(70) 
(71) 



The special value of /ii = 1 corresponds to the case where the potentials are exponential. We can now make 
some immediate progress: using (68) and (69) at the fixed point we show that 



a = /3;(1 - Hi) 



(72) 
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By further imposing ( 45 1 to ( 48 ) , we see that the powers in the potentials ( 70 ) and ( 71 1 are given by 

1 " ' 2 

3 



1 






1 


1 


- MP 




1 


1 


- Ms 




1 



1 - [l r 



R, 

Pj_ _ 


A 


a 


a 


Pp _ 


6 


a 


Q 


$9 = 


2 


a 


a 


Pr _ 


4 


a 


Q 



'— = - - 1 



We can also integrate ( 65 ) to give 

A* 



1 



1 — fj.i Hi — 1 



v 0i JA. 



« =3,P,9,r 



(73) 
(74) 
(75) 
(76) 

(77) 



1 - A*i 

A. a = fixed point: vanishing spatial curvature 

We shall begin by studying the fixed point a = because it corresponds to the situation where spatial 
curvature is sub-dominant, and so we expect the early behaviour to be closely matched by such fixed points 
even though a may not strictly vanish. Now (once again under the assumption that all the fii are well-defined 
and constant), we can make use of (77) and (68) to obtain 



V, 



Voii-ah/x)-^. 



(78) 
(79) 



We therefore have in terms of h and x once we have chosen a specific a, from which the fit are derived 
and hence the potentials determined. It allows us to replace the yi terms in the scale factor, and the scalar 
equations of motion at the fixed-point (where everything is constant), which leads to two equations for the 
variables, x and h. The Hamiltonian is not an independent equation but could have been used instead of 
say the scale factor equation of motion. 



Once we reach the fixed point, and using (45)- (48 1, and (69), the scalar field equation of motion (59) 
becomes 



(80) 
081) 



2xyj(Z + h) - 3ir 2 2/p(3 + h) - 2X g y g (2 + h) - 8y r (l + h) = 0, 
which in turn becomes upon using ( |79[ ) 

- 3(3 + h)V 0p {-ah) 2 - 2 ' a x 2 l a - 2V Qj ah{?, + h) - 8Vbr(l +h) + 4/i(2 + h)V og (-ah)- 1+2 / a x- 2/a 
We could equivalently write this as 

A(x 2/a ) 2 + B(x 2/a )+C = 0, (82) 

where 



A = -3{3 + h)V 0p (-ah) 2 - 2 / a , 

B = -2V 0j ah{3 + h)-8V 0r (l + h), 

C = 4h(2 + h)Vo g (-ah)~ 1+2 / a , 



(83) 
(84) 
(85) 



and solve the quadratic to find x(a, h, Voj/Vo p , Vo g /Vo p , Vor /Voy )- We also have the following for the scale- 
factor equation of motion ( 60 1 and Hamiltonian constraint ( 54 1 respectively 



(2h - 3)x 2 y 3 -.+ 3(2 - h)x 3 y p + 2(3 - 2h + 2h 2 )y g + 8(2 - K)xy r = p A , 
9x 2 yj - 15x 3 y p - 6y g [1 + X g x] - 2Axy r = -p A . 



(86) 
(87) 
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In particular the Hamiltonian ( 87 ) becomes a cubic equation in x 2 / a 

-A/a 



15a 3 h 3 Vc 



0/j 



ah 

x 



-6/c 



3ha[3haV 0j + 8V 0r 



ah 

x 



6(1 - 2h)V L 



ah 

x 



-2/c 



PA 



Now, any self tuning solution should be such that the parameters in the potential, a, Vo%, (i = j,p,r,g) are 
independent of pa, in other wor ds they should not be fine tuned for a particular value of p\. This implies 
that the solution for x 2 l a in (82 1 will be independent of p\ as there is no p\ term in it. However, this would 
then be inconsistent when substituted into the Hamiltonian constraint (88 1 where a p\ term is present. In 



order to avoid this potential inconsistency we must have that each of A, B and C va nish in (82) because it 

x 2/c 



and so x is allowed to depend on pa in ( 88 ) 



then does not give an equation for 

We can recover various results from ( 82 1 corresponding to different scale- factor evolutions (determined by 
h (38)) by considering the vanishing of various terms. We list some interesting examples in table III The first 



TABLE I: Examples of interesting cosmological behaviour for various fixed points with a = 0. 



Case 


cosmological behaviour 


VM 


v P (4>) 


VM 


V r (<t>) 


Stiff fluid 


H 2 oc l/o 6 












Radiation 


H 2 oc 1/a 4 


Cl</>° 









Curvature 


H 2 oc 1/a 2 













Arbitrary 


H 2 oc a 2h , h^O 


ci(l + h)cj>^- 2 








-^h(3 + h) Cl <j>^ 



three of these are consistent with the results of the previous section. However, perhaps the most interesting 
solution is that last one, labelled "Arbitrary". This corresponds to any cosmological evolution consistent 
with a power law expansion a oc including a matter dominated universe (h = —3/2), and inflationary 

expansion (—1 < h < 0). 

We will discuss the behaviour of these solutions in more detail in section [V] In particular we will switch 
the spatial curvature back on, and use numerical simulations to demonstrate how these various cosmological 
behaviours dominate at early times before giving in to self tuning at late times 2 , and so an asymptotically 
Milne Universe. We will also show that the most interesting of these solutions - the matter-like cosmology 
- suffers from a rapid gradient instability at the level of cosmological perturbations. This is unfortunate, 
but all is not lost. In the following section we obtain a new class of matter-like scaling solutions that exist 
for pA — 0. This may seem like a strange thing to do given that the Fab-Four was introduced to deal with 
a large cosmological constant, however, the analysis of the forthcoming section is simply a means to an end. 
As we will see in section [Vj in certain cases, these solutions will not dramatically alter their behaviour when 
Pa is turned on, and, more importantly, they are stable. 



B. a = fixed point with pa = 0: vanishing curvature and cosmological constant 



Setting pa = in the Hamiltonian constraint (11) implies that for fixed point solutions to exist the 
individual terms ( 12 )-( 15 ) no longer have to be time- independent. They can all scale with the Hubble 
parameter, for instance, as long as they scale in the same way. This then has a knock-on effect on the 
allowed solutions as we shall now see. In terms of the scaling parameters, the individual components of the 



2 Such late time behaviour only applies to matter and radiation, but not inflation 
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Hamiltonian (50)- (53 1 are replaced by 

Hp = -3x 3 y p (5~3cr 2 )H 6 - 3a -^, 

H g = -6y g [(1 - a 2 )H 2 -^ + X g xH 2 ~P g ~ a ~ lg } , 

H r = -24:xy r (l-a 2 )H 4 - a -P r , 



(89) 
(90) 
(91) 
(92) 



where the explicit H dependence is included. Now if we demand that all the individual terms in the 
Hamiltonian scale as H n (which, for constant h, is equivalent to H n ~ a nh ) then we obtain the following 
constraints, replacing (45)-(49) 



fa = 4 -2a- n, 
P P 

Pr 



lj = lp=lg = lr 



6 — 3a — n, 
2 - n, 
4 — a — n, 
—a. 



(93) 
(94) 
(95) 
(96) 
(97) 



The scalar held equations of motion |55|)-(58 1 are replaced by 



Ef'/{3a i H a+n ) = Qxy 3 (l - a 2 ) + 2(a + n- l)hxy 3 (l - a 1 ) + 2[xy 3 \\ - a 2 )] 1 - X 3 x 2 y 3 {\ - a 2 ), (98) 



E^/(3a 3 H a+n ) 
E^/{3a 3 H a+n ) 



-9x 2 y p (l - a 2 ) - 3(a + n- l)hx 2 y p {\ - a 2 ) - 3[x 2 y p (l - a 2 )]' 
+X p x 3 y p (l - a 2 ), 

-6X g y g (l -a)-2(a + n- l)hX g y g (l -a)- 2[\ g y g {\ - a)]' 
+2ng\lvgx{\ - a) + 2\ g y g {l - a) 2 , 



E^/(3a 3 H a+n ) = -24y r (l + h) 



-a 2 (l-a) + -(l-a 3 ) 



8y r 



-a 2 (l-a) + -(l 



the scale factor equations (|61[)- ( 64 ) become 



E) a> /(3a 2 H n ) = {-4x 2 y 3 + 2(1 - n)hx 2 y 3 - 2{x 2 y 3 )' + x 2 yj {l- a 2 )} , 

E^/{3a 2 H n ) = {2x 3 y p (3-<T 2 )-(l-n)hx 3 y p (3-a 2 ) + [x 3 y p (3- f j 2 )}'} 7 

E^/(3a 2 H n ) = {Ay' g + 4y g (2 - (1 - n)h) + 2(X g xy g )' - 2(1 - n)hX g xy g - 2y g (l + a 2 )} , 

Ei a) /{3a 2 H n ) = {l6xy r (l -a 2 )- 8(1 - n)hxy r (l - a 2 ) + 8[xy r {l - a 2 )]'} , 



(99) 
(100) 

(101) 



(102) 
(103) 
(104) 
(105) 



and, for completeness, the Hamiltonian constraint remains as before (recall we have ensured the time de- 
pendence in H factors out) 



3x 2 y 3 {3 - a 2 ) - 3x 3 y p (5 - 3a 2 ) - 6y g [(1 - a 2 ) + X g x] - 2Axy r (l - a 2 ) = 0. 



(106) 



The results obtained previously for the A's, y's and a (65)-(69) still apply, although of course with the 
new values for the P coefficients. We now look for a = fixed point solutions in order to understand the 
early-time behaviour. The Hamiltonian constraint ( 106 ) becomes 



9x 2 y 3 - I5x 3 y p - 6y g [1 + X g x] , -24xy r = 0, 



and the scalar equation of motion (98)- (101) becomes 



6xyj + 2(a + n — \)hxy 3 — X 3 x 2 y 3 — 9x 2 y p — 3(a + n — \)hx 2 y p 
+X p x 3 y p - 4X g y g - 2(q + n - l)hX g y g + 2[i g X 2 g y g x - 8y r (l + h) = 0. 



(107) 



(108) 
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We can follow the route taken in the previous subsection and use the fact that for all four potentials 

(109) 



(i = j,p,g,r) (68) and (69) still hold, giving 



a/Pi 



1 - Mi, 



which leaves ( 108 1 as 

[6 + (n + 2)h)]xyj - [9 + (3 + 2n)h]x 2 y p - 2X g y g [2 + h}- 8y r {l +h)=0. 



Now, if we take u = const then we can use (79 1 to rewrite the Hamiltonian (107) as 

4/ct 



6[1 + (n - 2)h]Vi 



(>;/ 



ah 

x 



[9a 2 h 2 V 0j + 24ahV 0r ] 



ah 



2/c 



15a 3 h 3 V 0p = 0. 



(110) 



(111) 



Similarly the scalar equation of motion (110) may be written as 

A/a 



2(2 + h)(2-n)hVo g (-^j ' + {[6 + (n + 2)h)}a 2 h 2 V aj ■+ 8(1 + h)ahV 0r } (—^ ' ( 112 ) 



+ [9 + (3 + 2n)h]a 3 h 3 V 0p = 0. 
If we now define the variables 



V j = a 2 h 2 V 0l , V 0p = a 3 h 3 V 0pi V 0g = V 0g , V 0r = ahV 0r , 



X 



{-ah/x) 2/a , 



then the Hamiltonian equation (111) and the scalar equation of motion (112) become 

- 6[1 + (n - 2)h]V 0g X 2 + [9V 0j + 24V 0r }X + 15% 
2(2 + h){2 - n)h% g X 2 + {[6 + (n+ 2)h]V 0j + 8(1 + h)V 0r } X + [9 + (3 + 2n)h]V 0p 



(113) 
(114) 

(115) 
(116) 



Our aim now is to find solutions of these equations. On the surface it would seem that we can just solve 
for X and n, and pick h and the Vo# at will. In particular, for a matter dominated epoch we would pick 
h = —3/2. However, it is worth pausing for a moment to think about our ultimate goal. We would like any 
such solution to be robust against switching /?a back on. To see how this might be possible, it proves useful 
to introduce 



hn, 



(117) 



since then, for constant h, we have % ~ H n ~ a" . So, if we can find solutions for positive h, the correction 
to the Hamiltonian from switching p\ back on, p\a~ n , may be expected to become less significant with 
expansion, meaning that the p\ — solution was indeed worth finding. At this point it also proves convenient 
to fix the freedom we have in re-defining the scalar field, by choosing the parameter a such that 



ah 



-1. 



We may then determine the powers of 4> that appear in the scalar potentials (70), (71 ) 

2 = h- (4/i + 2), 

3 = h -(6/1 + 3), 
■-n-2h, 

l = h- (4/i + l). 





1 


= Pj/<* = 


4 


— n 


1 


-to 




a 




1 


= P P /a = 


6 


— 11 


1 


- m p 




a 




1 


= P g /a = 


2 


— n 


1 


-Ms 




a 




1 


= p r /a = 


4 


— n 


1 


- fl r 




a 



(118) 

(119) 
(120) 
(121) 
(122) 
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We shall now focus on reproducing a matter epoch, so we plug h 
This gives 



-3/2 and a = 2/3 into ([Tl5| and (fUgf 



6(n + 4)V 09 A 2 + [9F 0j -24Vb,.]A-15Uop = 

" 2 4- [Cn -4- .Wn, -4- 4t^_l X - ( - 



-{ft + 3)V 0g X 2 + [{h + 3)V 0j + AV Qr ] X-[- + 2n)V 0p = 



(123) 
(124) 



We want these two equations to have a common root. Although we will miss the full spectrum of possibilities, 
it is convenient and cleaner to impose the stronger constraint that the two equations have two common roots. 
This requires that the following equations hold for some constant \ 



6 X {n + 4)y 0s 

X [W 0j - 2W 0r ] 



{n + 3)V 0g , 

{h + 3)V 0j + W Qr , 

9 2n ] V 0p . 



(125) 
(126) 

(127) 



Now if Vq p , Vo g 7^ then it follows from ( 125 1 and (1271 that n = — |, — J. However, as explained above, we 



are only interested in h > 0, so we assume that one 



of 



^Opj Vo g vanishes. The two possibilities are presented 



in table [II] In each case, either Paul or George is switched on/off, and the interplay of John and Ringo gives 



TABLE II: Examples of matter-like cosmological fixed points with a = and pa = 0. 



Case 


cosmological behaviour 






v B (<t>) 


K(4>) 


Matter I 


H 2 oc 1/a 3 




C2<f> A+(> 





2n — 3 „ 
16(2n + 7)(n + 6) Cl< P 


Matter II 


H 2 oc 1/a 3 









(n + 3)(2n + 5) ,h + 6 
8(2n + 7)(n+0) Ll ^ 



rise to a cosmology that mimics a matter dominated epoch. This is true for any value of n, and in particular, 
for fi > 0. We will see in the next section that these solutions are indeed robust against switching on 
and there exists large regions of parameter space (n,ci,C2) for which their cosmological perturbations are 
stable. 



V. SUMMARY OF COSMOLOGICAL SOLUTIONS 



The goal of this paper was to establish whether or not the Fab-Four could, in principle, accommodate a 
consistent cosmological history. In particular, is the Fab-Four consistent with an early period of inflation, 
followed by a radiation and matter dominated epoch, during which nucleosythesis takes place and structures 
begin to form? At late times, we want another period of inflation, before self-tuning kicks in and the Universe 
enters a late time Milne solution, whatever the value of the cosmological constant. Our dynamical systems 
analysis, allied with the numerical solutions to be presented shortly, demonstrates that each desired epoch 
can be individually realised with a judicious choice of potentials. It is not too difficult to imagine that 
one could, in principle, combine the various choices in such a way that one particular choice dominates the 
dynamics at one particular epoch, thereby reproducing the desired cosmic history in its entirety. Recall that 
the cosmological constant is always assumed to be large, dominating over any other sources. Including extra 
sources explicitly is straightforward as we know the scaling behaviour of both radiation and matter with 
scale factor evolution. 

Having obtained desirable solutions, we should also demonstrate that they are stable to perturbations, or 
at least stable enough that they will survive for the e-foldings necessary for structures to form. Scalar-tensor 
theories can be plagued with ghost and gradient instabilities (see [4] for a review of the subject), which is 
one of the main reasons it is proving so challenging to develop successful modified theories of gravity. This 
involves lengthy calculations in each case, and for that reason we restrict attention, for now, to the epoch of 
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greatest interest - matter domination. We will see that the matter epoch is perturbatively stable for a large 
class of cases. 

We will now summarise how each desired epoch can be reproduced within the Fab-Four, reinforcing our 
analytic calculations with numerical simulations of the full system. In each case we will plot the evolution 
of the deceleration parameter, q — —ad /a 2 , which is expected to give q < (inflation), q — 1 (radiation), 
q = 0.5 (matter) and q = (self-tuning). As promised, the matter epoch will also include a discussion of 
stability. 



A. The inflationary epoch 



It is natural ask how any self-tuning scenario can accommodate inflation, be it early or late. The point is 
that self-tuning can be a late time attractor, and that inflation can happen beforehand. Clearly this would 
be perfectly consistent with what we observe. In section IV A we saw how power law inflation, a ~ t p , with 
p > 1 could be achieved for arbitrary p\ with the choice: 



W) = ci 



P-1 

P 



v p {<i>) = v g {4>) = o, 



3p — 1 

16" 



(128) 



where we have used the fact that p = —1/h, and have fixed the freedom to redefine the scalar field by setting 
ah = — 1. For vanishing spatial curvature, the field equations (22 1 and (231 can be straightforwardly solved 
to give 



<p = va, 



H 2 = 



PPA 



3 Cl (3p+l 



-{va) 



-2/p 



(129) 



where v is an arbitrary constant. The inflationary case actually presents a qualitatively different behaviour 



to the matter and radiation solution, as could have been guessed from (|67|). In (67) we see that one may 



expect a — to be a repeller fixed-point for 1 + h < 0, such as matter and radiation, but an attractor for 
1 + h > 0. Although this is rather naive, as there are many more variables to consider, this is indeed what 
is observed. Switching on spatial curvature and performing the full evolution drives the system to the a = 
inflationary solution, as seen in Fig. [T] Of course, in a more complicated scenario in which potentials are 



1r 



0.5 



-0.5 




10" 



10" 



a(t)/a(t f ) 



10" 



10 u 



FIG. 1: Plot of the deceleration parameter, q, (lower curve) and 100a (upper curve) for the inflationary 
case for p = 3. The parameters used were k = —10, pa = 1000, (^initial = 0.1, c\ = 1.5 



"sewn together" so that inflation gives way to radiation one expects self-tuning to dominate at very late 
times. 



15 



B. The radiation epoch 



In section IV A we found that for arbitrary p\ the following potentials will mimic a radiation dominated 
Universe, 



y J (0) = Cl 6 , v p (<f>) = Q, 

with a vanishing spatial curvature solution of 

<b = z/a, H 2 = 



Cl 



1± Jl 



15c 2 



(va)- 



(130) 



(131) 



where v is an arbitrary constant. Note that as c\ — > the negative root above has a well defined limit. This 
corresponds to the case where only George is turned on. We could easily have guessed that this would mimic 
radiation since the scalar equation of motion imposes the constraint R — 0. This is equivalent to saying that 
the trace of the effective energy momentum tensor vanishes. 

When we reintroduce the spatial curvature, self-tuning kicks in. This is explicitly demonstrated by the 
numerical solutions presented in figure [2j 

1.4i 1 

1.2 

1 



0.8 
q(t) 

0.6 
0.4 
0.2 

i 



10- 4 10- 2 

a(t)/a(t f ) 



10" 



FIG. 2: Plot of the deceleration parameter for an initially radiation-like solution, with k = —0.1, 

PA = 1000, ^initial = 0.1, Ci = -0.5. 



C. The matter epoch 



We finally turn to the matter epoch. We have identified three classes of potential that can mimic a matter 
dominated universe. Let us consider the first two of these, dubbed "Matter I" and "Matter II" , and identified 
in section [IV B| for vanishing curvature and cosmological constant. The "Matter I" potentials have vanishing 
V g , and are given by 



v j {4>) = Cl( j> 



n+4 



v p {4>) = c 2 <f> 



71 + 6 



v g (<p) = o, 



Vrift>) 



2h - 3 



16(2n + 7)(n + 6) 

For vanishing spatial curvature and cosmological constant, this admits a solution, 



ci 



,n+6 



(f) = va, 



H 2 = 



2h 



2{2n + 7) 



ci 

C2 



{va) 



-3 



(132) 



(133) 



16 



where v is an arbitrary constant. The "Matter II" potentials have vanishing V p , and are given by 

V^) = C X 6^+\ V p (6) = 0, V g (<t>) = C 2 <t> A+3 , V r {4>) 



(n + 3)(2n + 5) fl+6 



8(2n + 7)(n + 6) 

This time, for vanishing spatial curvature and cosmological constant, we have the solution, 



H 



2(2h + 7) 



2n + 9 



C2 
Cl 



(134) 



(135) 



where ^ is an arbitrary constant. Note that this solution is not well defined as c\ — > 0, when we expect to 
recover a radiation-like Universe. 

Of course, these solutions are only valid when = 0, which goes against the spirit of the Fab-Four and 
self-tuning. What happens when we switch p\ back on? To get an idea, let us perturb about these solutions 
and compute the Hamiltonian (22 1 and scalar equation of motion (23) to leading order. To this end, we set 



6 = ua + 86(N), 



H = 



8H 2 (N), 



with 



f 


2n+9 




2(2n+7) 




"2(2*1+7) 




2n+9 



Matter II 



and expand (22) and (23) to leading order in 8. The result is 

SUj + SH P + 5H g + SH r = -pa- 



+ sty*) + sifyti + SEW = 0. 



where 



JH 2 ,„ „.<ty n 
2^ + ^ + 6)^ + 2 

SH 2 , . 86 
3^ + (, + 9)^ + 3 



86\' 



5H 



fin 



rtr 



5H 2 .86 

8H 2 
2 1P 



n + 3 
n + 4 



+ (n + 6)-f- + 



(136) 



(137) 



(138) 
(139) 



(140) 
(141) 
(142) 
(143) 
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and 



SE 



(e^^/Hc 



Ha ] = 



3 



. SH 2 \ , JH 2 
3 1 — 1 + 2(n + 3)- 



H 2 



+ (n + 3)(n + 6)- 
15 f5H 2 



H 2 



+ 3(n + 4) 



(144) 



H 2 
9 



3 n 



9\ SH 2 



H 2 



An 



63 

T 



6EM =5[Ef>d- ff" 



SE^ =S(E^4>/Hc 



Hg 

2 

~2 



n 



n + A 
SH 2 



(ft + 9) 



H 2 



+ (n + 3)- 



ff 2 



(n + 6)- 



(145) 
(146) 

(147) 



Note that the unperturbed Hamiltonian densities vary as H oc a™, which suggests that even a large cosmo- 
logical constant may be regarded as sm all p erturbation, provided p\/d n <C 1. We now use (138) to find 
5H 2 /H 2 , and substitute the result into (139). In each case we obtain 



= f(n;c 1 ,c 2 ) 



where 



/(n;ci,c 2 ) 



This is straightforwardly solved to give 



2(2n+27))(2w+7) 3 C2 
' (2n+9) 2 (4n 2 +28fi+9)c^ 

(2ra+9) 2 (2n 2 +5n-9)ci 
' 16(2n+7)(4ri 3 +32n 2 +71n+27)c 



PA 

Matter I 
Matter II 



A, 



A? 



2/ PA 



a" 1 2 9n [va) n 

where A\ and A 2 are integration constants. It immediately follows that 

SH 2 „ B 2 



i? 2 



Si 



PA 

(i/o)"' 



(148) 



(149) 



(150) 



(151) 



where the constants Bi,B 2 and g are related to Ai,A 2 and / respectively. It is now clear that for n > 0, 
Ci , C2 7^ 0, the solutions ( 133 1 and ( 135 ) are late time attractors for vanishing curvature, even when we switch 
on pa- This view is reinforced by the plots shown in figures [3] and [4] In each case, the evolution mimics 
a matter dominated epoch for a long time even in the presence of a cosmological constant. When we also 
include curvature in the numerical simulation, the solution ultimately gives in to self tuning and asymptotes 
to a Milne universe. As we will see in the next section, one can choose h > 0, C\ and c 2 such that we have 
stable cosmological perturbations for Matter I and Matter II. 

The third class of potentials to mimic matter domination, which we dub "Matter III" , are given by 



v p {<t>) = W) = o, v r {4>) 



16 



(152) 



and follow from the " arbitrary" row of table |I| giving a matter-like solution for arbitrary p\ when h = —3/2 
and a = 0. The explicit solution is 



2 Pa 
81ci 



(153) 



18 




FIG. 3: Plot of the deceleration parameter for the Matter I solution, with k = —10 , 4>i n itiai — 3.0, 
Pa = 1, h — 6, ci = 1, C2 = 10, consistent with stable cosmological perturbations, as shown in section VC 1 




10 
a(t)/a(t f ) 



FIG. 4: Plot of the deceleration parameter for the Matter II solution, with k = —0.1, 



iai — 0.1, pa — 1, 



h = 2, c\ = 1, C2 — 8, consistent with stable cosmological perturbations, as shown in section VC 1 



where v is an arbitrary constant. The numerical solution with non- vanishing curvature is shown in figure 
[5j Whilst this solution seems well behaved at the level of FLRW cosmology, we will see that it is unstable 
under cosmological perturbations. 

In summary, then, we have found three classes of potential that mimic a matter dominated cosmology, 
even when the source is dominated by a cosmological constant. For the Matter I and Matter II potentials, we 
can choose parameters such that the solutions are stable under cosmological perturbations. These solutions 
are remarkable. The Fab-Four potentials have forced the scalar to screen the pressure component of the 
cosmological constant before it screens the energy density. This allows for an intermediate, and in some 
cases, pathology-free, period resembling a matter dominated cosmology, even in the absence of a pressureless 
source. At the latest times, the energy component of the cosmological constant is also screened and the 
solutions evolve to an asymptotically Milne universe, which is, of course, equivalent to a patch of Minkowski 
space. 
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FIG. 5: Plot of the deceleration parameter for an initially Matter III solution, with k = —0.1, p\ — 1000, 

^initial = 0.1, C\ = 1. 



Stability during matter epoch 



As promised, we will now consider cosmological perturbations about each class of matter solution. Since 
we are interested in the phase prior to self-tuning, we shall restrict attention to vanishing curvature, and 
the solutions given by (133), (135) and (153). Vacuum perturbations about spatially flat cosmologies were 
studied in detail for Horndeski's theory in [10 . Working in unitary gauge, the scalar takes on its background 
value, 4> = 4>(t), whereas the line element is given by 



ds 2 = -~N 2 dt 2 + ^{dx 1 + N l dt){dx j + N j dt), 

with 

iV = l + d, Ni = dip, 7tf =a 2 (t)e 2C (5 ij +h ij ). 
The quadratic action for tensor perturbations is found to be 

1 



St, 



dtd i xa A 



Tr 
a 



a (VM 2 



(154) 
(155) 

(156) 



where Gt and -Ft depend on the potentials evaluated on the background solution. Note that [TD] uses the 
DGSZ form of Horndeski's theory [5], so Gt and Tt are given in terms of DGSZ potentials. The DGSZ 
potentials for the Fab-Four are given in Appendix C of [T3], so we can use this to extract the form of Gt 
and Tt for the cases we are interested in. We shall spare the reader the details, since they add little to this 
discussion. 

Upon integrating out the lapse and the shift, the scalar perturbations yield the following effective action, 



S. 



scalar 



= / dtd A 



GsC 2 



T s 



(VC) 1 



(157) 



Again, the coefficients C?s and Ts depend on the DGSZ potentials evaluated on the background solution. 
Now, to avoid a ghost instability, we require 



Gt > 0, Gs > 0, 
and to avoid a gradient instability, we require 

T T > 0, T s > 0. 



(158) 
(159) 
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We first consider the Matter I and Matter II solutions given by (1331 and (135), which, as we have seen, 
are late time cosmological attractors for vanishing spatial curvature, even when the cosmological constant is 
non-zero. Plugging this solution into the appropriate formulae for Qx, etc, we find that 



Gt,s oc Cit 3 



in+2 



T T ,s oc cii5 fl+2 . (160) 

However, the crucial thing to know is the sign of each of these coefficients. These are plotted in figure [6] In 
each case the wrzshaded area correspond to regions of parameter space in which all coefficients are positive. 
We immediately see that it is not too difficult to guarantee the absence of ghost and gradient instabilities 
for both Matter I and Matter II. 
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FIG. 6: Stable and unstable regions of parameter space for Matter I solution (left panel) and Matter II 
solution (right panel). In each the shaded regions have one or more ghost/gradient instability. The white, 
unshaded, regions have no such pathologies and are stable against cosmological perturbations. 



The Matte r III potentials do not give rise to such well behaved cosmological perturbations about the 
solution (153). Tensor perturbations are fine, with no ghost or gradient instability 



(161) 



Q T - 28M 4 t 2 , T T - 104M 4 i 2 , 
for some real mass scale M. In contrast, the scalar perturbations exhibit a gradient instability 

S s ~^M 4 i 2 , F S ~-~MH*. (162) 

For modes of wavelength A this instability manifests itself on timescales ti nsta biiity ~ -Vl c s| ; where \c s \ — 
^J- ~ 0(1). This is far too quick, and rules out the Matter III solution as part of a viable cosmology. 
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VI. CONCLUSIONS 



Obtaining a sensible cosmology in the presence of a large and changing contribution to the vacuum energy 
Pa is one of major challenges facing the self tuning scenario that we have developed in [21 E] . A standard 
cosmological evolution arising out of General Relativity with large p\ would be totally unacceptable apart 
from perhaps in the early Universe where it would drive a period of accelerated expansion. There would be 
no way of exiting this period of inflation and obtaining a period of radiation and matter domination required 
for nucleosynthesis and structure formation. The goal of this paper is to demonstrate that for the Fab-Four , 
it is indeed possible to obtain a sensible cosmological history even in the presence of a large p\ contribution 
at all times. In other words, the self tuning of the cosmological constant can be accommodated in a sensible 
cosmological timeline. 

To show this we have developed a dynamical systems approach in which fixed point solutions corresponding 
to inflationary, radiation and matter dominated solutions are made manifest. Two key approaches are 
developed. In the first, we explicitly include the p\ contribution and by demanding that all contributions 
in the Hamiltonian constraint remain independent of the Hubble parameter, we show that there exist a 
class of scaling solutions corresponding to the cosmologies we are looking for (recall there is no requirement 
here to actually include matter or radiation sources, the scalar field is doing the work for us). However, it 
turns out that this particular matter dominated solution (called "Matter III"), whilst perfectly acceptable 
at the background level, actually contains a gradient instability when perturbed, an instability that would 
grow on too fast a timescale to be compatible with observations. This has led us to consider a second 
complementary approach. Rather than obtain background solutions in the presence of p\ we set it to zero 
and look for consistent solutions that can also mimic matter domination (H 2 cx a -3 ). Our requirement that 
the Hamiltonian constraint be independent of the Hubble parameter is now lifted and this allows for more 
freedom introducing an extra parameter (we call n or equivalently, n) in the background solutions. We find 
two new classes of matter-like solutions (called "Matter I" and "Matter") for arbitrary n. 

It would seem to go against the self-tuning spirit of the Fab-Four that the "Matter I" and "Matter II" 
solutions only correspond to fixed points for vanishing cosmological constant. However, we have shown that 
they can still represent an excellent approximation even when a large p\ is turned on, provided n > 0. 
This is because the solution gets corrected by p\/a n , a correction that decreases with time as the scale 
factor grows. Using analytic methods, we also showed that for vanishing spatial curvature, these solutions 
are cosmological attractors for h > 0. Once spatial curvature is reintroduced alongside the cosmological 
constant, we are forced to use numerical simulations which reproduce the expected behaviour: a long period 
of matter-domination, before asymptoting to the self-tuning Milne Universe. For the record, the Lagrangians 
for the "Matter I" and "Matter II" cosmologies are given by 

2fi — 3 
16(2n + 7){n + 6) 

£"Matt cr ii" = c^+^v^v^ + c 2 r +3 R - Cl r +6 g, 

8{2n + 7){n + 6) 

with the parameters n > 0, C\ and c 2 chosen so that they correspond to the perturbatively stable unshaded 
regions shown in figure [6] Indeed, perhaps the most exciting new feature of the corresponding solutions, is 
that they can be stable against cosmological perturbations for large regions of parameter space when spatial 
curvature is negelcted. This opens up the possibility of a sensible matter dominated period of evolution, 
hence of structure formation in the Fab-Four scenario. Furthermore, these solutions are behaving in such a 
way that the scalar screens the pressure component of the cosmological constant before the energy density. 
At least for homogeneous and isotropic backgrounds, this suggests that the cosmological constant is being 
forced to behave like cold dark matter. It is certainly tempting to ask whether such behaviour extends to 
inhomogencous solutions, and recent results suggest that it may well be possible to have a Fab-Four scenario 
satisfying current solar system constraints [18] . 

There is much that remains to be done. We have not yet obtained a full cosmology, but the fact that we 
have a class of background polynomial potentials that we know can provide the various cosmological epochs 
we want to reproduce offers us some direction. Indeed we can speculate as to how we might sew together these 
interesting potentials to achieve the desired results. The point is that the scalar field is continually evolving, 
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so we could arrange for the potential to correspond to different fluid behaviours for different ranges of <fi. 
For example, if we want radiation domination for H 2 > H 2 q and "Matter I" like behaviour for H 2 < H 2 q , 
we might propose a Lagrangian of the form 



4 

£ = TpT^R + "Matter I" 



ci-0(l),C2~0(m- 2 ) 



(163) 



Plugging in the "Matter F solution (133), we see that 



'"' , .*R~m 4 O fl (^ , £"MatterP ~ m V\ (164) 



H 2 



where a eq = (m/ H eq ) 2 / 3 is the value of the scale factor when H = H eq . We see that the "Matter I" terms 
dominate for a > a eqi as desired for mater domination. For a < a eq the (j) n R term becomes important, and 
might be expected to dominate the dynamics, yielding an earlier period of radiation domination. The same 
may be done to evade the graceful exit problem due to the inflationary solutions being attractors. 
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